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Abstract 

A remarkable fact about spherically-symmetric neutron stars in hydrostatic equilibrium — the 
so-called Schwarzschild stars — is that the only physics that they are sensitive to is the equation 
of state of neutron-rich matter. As such, neutron stars provide a myriad of observables that may 
be used to constrain poorly known aspects of the nuclear interaction under extreme conditions 
of density. After discussing many of the fascinating phases encountered in neutron stars, I will 
address how powerful theoretical, experimental, and observational constraints may be used to place 
stringent limits on the equation of state of neutron-rich matter. 
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I. INTRODUCTION 



A neutron star is a gold mine for the study of the phase diagram of cold baryonic matter. 
While the most common perception of a neutron star is that of a uniform assembly of 
neutrons packed to densities that may exceed that of normal nuclei by up to an order of 
magnitude, the reality is far different and significantly more interesting. Indeed, the mere 
fact that hydrostatic equilibrium must be maintained throughout the neutron star, demands 
a negative pressure gradient at each point in the star; otherwise the star would collapse under 
its own weight. This model- independent fact yields nuclear densities — at least for most 
realistic equations of states — that span over 11 orders of magnitude; from 10 4 to 10 15 g/ cm 3 . 
Recall that in this units nuclear-matter saturation density equals po — 2.48 x 10 14 g/cm 3 . 
What novels phases of baryonic matter emerge under these conditions is both fascinating 
and unknown. Moreover, most of the exotic phases predicted to exist in neutron stars can 
not be realized under normal laboratory conditions. Whereas most of these phases have a 
fleeting existence here on Earth, they become stable in neutron stars as a consequence of 
the presence of enormous gravitational fields. 

To establish the fundamental role played by the equation of state on the structure of 
spherically-symmetric neutron stars in hydrostatic equilibrium, we start with the Tolman- 
Oppenheimer-Volkoff (TOV) equations — the extension of Newton's laws to the domain 
of general relativity. The TOV equations may be expressed as a coupled set of first-order 
differential equations of the following form: 
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where G is Newton's gravitational constant, while P(r), £(r), and M(r) represent the 
pressure, energy density, and enclosed-mass profiles of the star, respectively. Note that the 



last three terms (enclosed in square brackets) in Eq. (la) have a general- relativistic origin 



Remarkably, the only input that neutron stars are sensitive to is the equation of state of 
neutron-rich matter. Indeed, changes in pressure and enclosed mass as a function of radius 
(left-hand side of the equations) depend not only on the values of these quantities at r, but 
also on the "unknown" energy density £(r) of the system. Thus, no solution of the TOV 
equations is possible until an equation of state (i.e., a P vs £ relation) is supplied. 

In this manuscript we discuss the various fascinating phases of baryonic matter that are 
predicted to exist in neutron stars, but inaccessible under normal laboratory conditions. 
After briefly discussion the theoretical models used in this contribution, we focus on recent 
theoretical, experimental, and observational constrains that place stringent limits on the 
equation of state of neutron-rich matter. 



II. ANATOMY OF A NEUTRON STAR 

Neutron stars contain a non-uniform crust above a uniform liquid mantle. See Fig. [I] for 
what is believed to be an accurate rendition of the structure of a neutron star. 
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A. The Outer Crust 



The outer crust is understood as the region of the star spanning about 7 orders of magni- 
tude in density; from about 10 4 g/cm 3 to 4 x 10 11 g/cm 3 pQ. At these densities, the electrons 

- which are an essential component of the star in order to maintain charge neutrality - 
have been pressure ionized and move freely throughout the crust. Moreover, at these "low" 
densities, 56 Fe nuclei arrange themselves in a crystalline lattice in order to minimize their 
overall Coulomb repulsion. This is the structure of the outermost layer of the crust. How- 
ever, as the density increases (and one moves away from the surface of the star) 56 Fe is no 
longer the most energetically favorable nucleus. This is because the electronic contribution 
to the energy increases faster with density than the nuclear contribution. As a result, it 
becomes energetically advantageous for the energetic electrons to capture on the protons 
and for the excess energy to be carried away by neutrinos. The resulting nuclear lattice is 
now made of nuclei having a neutron excess larger than that of 56 Fe. As the density con- 
tinues to increase, the nuclear system evolves into a Coulomb lattice of progressively more 
neutron-rich nuclei until a "critical" density of about 4 x 10 11 g/cm 3 is reached. At this point 
the nuclei are unable to hold any more neutrons; the neutron drip line has been reached. 




FIG. 1: A rendition of the structure and phases of a neutron star (courtesy of Dany Page). 



B. The Inner Crust 



The inner crust of the neutron star comprises the region from neutron-drip density up 
to the density at which uniformity in the system is restored (approximately 1/3 to 1/2 of 
normal nuclear matter saturation density). At these densities the system exhibits rich and 
complex structures that emerge from a dynamical competition between short-range nuclear 
attraction and long-range Coulomb repulsion. At the lower densities present in the outer 
core, these length scales are well separated and the system organizes itself into a crystalline 
lattice of neutron-rich nuclei. In contrast, at a much higher density of the order of half 
of nuclear-matter saturation density, uniformity in the system is restored and the system 
behaves as a uniform Fermi liquid. Yet the transition region from the highly-ordered crystal 
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to the uniform liquid mantle is complex and not well understood. Length scales that were 
well separated in both the crystalline and uniform phases are now comparable, giving rise to 
a universal phenomenon known as "Coulomb frustration". It has been speculated that the 
transition to the uniform phase must go through a series of changes in the dimensionality and 
topology of these complex structures, colloquially known as "nuclear pasta" [21 E]- In Fig.|2]a 
snapshot obtained from Monte- Carlo/Molecular-Dynamics simulations of a nuclear system 
at densities relevant to the inner crust are displayed (H E]. The figure displays how the 
system organizes itself into neutron-rich clusters (i.e., "nuclei") of complex topologies that 
are surrounded by a vapor of (perhaps superfluid) neutrons. Such complex pasta structures 
may have a significant impact on various transport properties, such as neutrino and electron 
propagation. 




FIG. 2: (color online) A snapshot of a Monte Carlo simulation for a configuration of 4,000 nucleons 
at abaryon density of n = 0.025 fm~ 3 (a sixth of normal nuclear matter saturation density), a proton 
fraction of Y P = Z/A = 0.2, and a temperature of T=l MeV. 



C. The Stellar Core 

As the density continues to increase, the neutron-rich nuclei will "melt" and uniformity 
in the system will be restored. At these densities (of the order of l/3-to-l/2 of normal 
nuclear-matter saturation density) the naive perception of a neutron star, namely, a uni- 
form assembly of closely-packed neutrons will be realized, albeit for the presence of a small 
percentage (of the order of 10%) of protons and electrons required to maintain chemical 
equilibrium and charge neutrality. Although the non-uniform crust displays fascinating and 
intriguing dynamics, its structural impact on the star is rather modest. Indeed, the crust 
accounts for about 10% of the radius of the neutron star and for only a few percent of its 
mass. Most of the mass of the neutron star is contained in its uniform core. Yet the high- 
est density attained in the core depends critically on the equation of state of neutron-rich 
matter which at those high densities is poorly constrained. The cleanest constraint on the 
equation of state at high-density will emerge as we answer one of the central questions in 
nuclear astrophysics: what is the maximum mass of a neutron star? Or equivalently, what 
is the minimum mass of a black hole? Note that if the equation of state is "soft" , very high 



4 



densities may be reached in the stellar core. At such high densities new states of matter may 
develop as the quarks within the hadrons become deconfmed. Such an exciting possibility 
will not be considered further in this manuscript. 



III. CONSTRAINTS ON THE EQUATION OF STATE 



Before addressing the role that recent observables play in constraining various theoretical 
description of the equation of state, we introduce the relativistic mean-field models that are 
used to compute these observables. 

Relativistic mean-field descriptions of the ground-state properties of medium to heavy 
nuclei have enjoyed enormous success. These highly economical descriptions encode a great 
amount of physics in a handful of model parameters that are calibrated to a few ground-state 
properties of a representative set of medium to heavy nuclei. An example of such a successful 
paradigm is the relativistic NL3 parameter set of Lalazissis, Ring, and collaborators [HI [7] . 

The Lagrangian density employed in this work is rooted on the seminal work of Walecka, 
Serot, and their many collaborators (see Refs. [HI 13 [10] and references therein). Since 
first published by Walecka more than three decades ago [H], several refinements have been 
implemented to improve the quantitative standing of the model. In the present work we 
employ an interacting Lagrangian density of the following form [TTJ [T21 H3] : 



^int = i> 



2 p 2 

- ^(gs^-^r+^gXv 1 ") "+A v (^b, • b") (fw) ■ (2) 

The original Lagrangian density of Walecka consisted of an isodoublet nucleon field (ip) 
together with neutral scalar (0) and vector (V 11 ) fields coupled to the scalar density ('ipip) 
and conserved nucleon current {ip^ip) , respectively [H|. In spite of its simplicity (indeed, 
the model contains only two dimensionless coupling constants), symmetric nuclear matter 
saturates even when the model was solved at the mean- field level [8j. By adding addi- 
tional contributions from a single isovector meson (6 M ) and the photon (A^), Horowitz and 
Serot [H] obtained results for the ground-state properties of finite nuclei that rivaled some 
of the most sophisticated non-relativistic calculations of the time. However, whereas the 
two dimensionless parameters in the original Walecka model could be adjusted to repro- 
duce the nuclear saturation point, the incompressibility coefficient (now a prediction of the 
model) was too large (K > 500 MeV) as compared with existing data on breathing-mode 
energies [13]. To overcome this problem, Boguta and Bodmer introduced cubic (k) and 
quartic (A) scalar meson self-interactions that accounted for a significant softening of the 
equation of state (K = 150 ±50 MeV) [16]. Two parameters of the Lagrangian density of 
Eq. ([2]) remain to be discussed, namely, ( and A v . Both of these parameters are set to 
zero in the enormously successful NL3 model, suggesting that the experimental data used in 
the calibration procedure is insensitive to the physics encoded in these parameters. Indeed, 
Miiller and Serot found possible to build models with different values of ( that reproduce 
the same observed properties at normal nuclear densities, but which yield maximum neutron 
star masses that differ by almost one solar mass [12]. This result indicates that observa- 
tions of massive neutron stars — rather than laboratory experiments — may provide the 
only meaningful constraint on the high-density component of the equation of state. Finally, 
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the isoscalar-isovector coupling constant A v was added in Ref. [13J to modify the density 
dependence of the symmetry energy. It was found that models with different values of A v 
reproduce the same exact properties of symmetric nuclear matter, but yield vastly different 
values for the neutron skin thickness of heavy nuclei and for the radii of neutron stars [17J. 
The Parity Radius Experiment (PREx) at the Jefferson Laboratory promises to measure 
the skin thickness of 208 Pb accurately and model independently via parity-violating electron 
scattering [181 US]- PREx will provide a unique experimental constraint on the density de- 
pendence of the symmetry energy due its strong correlation to the neutron skin of heavy 
nuclei [20]. 

A. Theoretical Constraints 

One of the most stringent constraints on the equation of state of low density neutron- 
rich matter emerges from theoretical considerations, namely, from the universality of dilute 
Fermi gases with an "infinite" scattering length (a). In this limit the only energy scale in 
the problem is the Fermi energy (ep) 5 so the energy per particle is constrained to be that of 
the free Fermi gas up to a dimensionless universal constant (£) that is independent of the 
details of the two-body interaction [2T]. That is, 



FIG. 3: (color online) Equation of state of pure neutron matter as a function of the Fermi momen- 
tum. Predictions are shown for the accurately calibrated NL3 j6[ [7J (green line) and FSUGold 
(blue line) parameter sets. Shown also are various microscopic descriptions — including a model- 
independent result based on the physics of resonant Fermi gases by Schwenk and Pethick [25] (red 
region) . 
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To date, the best theoretical estimates place the value of the universal constant around 
^ 0.4 [2H [2211231 [21]. 

Although the neutron-neutron scattering length is large indeed (a nn = —18.5 fm), pure 
neutron matter deviates from unitarity due to a non-negligible value of the effective range 
of the neutron-neutron interaction (r e = +2.7 fm). Thus, corrections to the low-density 
equation of state of pure neutron matter must be computed for kp^r' 1 ~ 0.4 fm -1 . Such 
corrections have been recently computed by Schwenk and Pethick [25] , with their results dis- 
played as the red hatched region in Fig [3} Also shown are the predictions of two microscopic 
models based on realistic two-body interactions, one of them being the venerated equation of 
state of Friedman and Pandharipande [26J. Finally, the predictions of NL3 and FSUGold are 
also shown. It is gratifying that the softening of the symmetry energy of FSUGold — caused 
by incorporating constraints from breathing-mode energies [TT] - appears consistent with 
the physics of resonant Fermi gases. Such a powerful universal constraint should be routinely 
and explicitly incorporated into future determinations of density functionals. Indeed, such 
a constrain appears to rule out many of the models displayed in Fig. 2 of Ref. [20J. 

B. Experimental Constraints 




FIG. 4: (color online) Binding energy per nucleon as a function of baryon density (expressed in units 
of the saturation density po = 0.148 fm -3 ) for symmetric nuclear matter. Theoretical predictions 
are shown for the NL3 |6l [7] (green line) and FSUGold [11] (blue line) models. Shown in the inset 
is a comparison between the equation of state extracted from energetic nuclear collisions |27J and 
the predictions of these two models. 

Energetic nuclear collisions may be used to constrain the high-density behavior of nucle- 
onic matter. To illustrate this point we display in Fig. [4] the binding energy per nucleon 
of symmetric nuclear matter as a function of the baryon density as predicted by both the 
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NL3 and FSUGold models. Note that both models reproduce the equilibrium properties of 
symmetric nuclear matter and display the same quantitative behavior at densities below the 
saturation point. Yet their high-density predictions are significantly different. This emerges 
from a combination of two factors. First, FSUGold predicts an incompressibility coefficient 
K considerably lower than NL3, namely, 230 MeV vs 271 MeV. Second, and more impor- 
tantly, FSUGold includes an omega-meson self-energy coupling [labeled by ( in Eq. (|2])] 
that is responsible for a significant softening at high density. We now compare the pre- 
dictions of these two models against results obtained from energetic nuclear collisions that 
can compress baryonic matter to densities as high as those predicted to exist in the core of 
neutron stars. The inset in Fig. [4] provides us with such a comparison. By analyzing the 
manner in which matter flows after the collision of two energetic gold nuclei, the equation of 
state of symmetric nuclear matter was extracted up to densities of 4-to-5 times saturation 
density [27] . Figure [I] seems to rule out overly stiff equations of state (such as NL3). And 
while it is gratifying that FSUGold is consistent with this analysis, one must stress that 
the connection between energetic nuclear collisions and the equation of state of cold nuclear 
matter is model dependent. 

C. Observational Constraints 

A recent observation that seems to suggest a hard equation of state is that of the low- 
mass X-ray binary EXO 0748-676. Note that such a binary system consists of a neutron 
star accreting mass from a normal (non-compact) companion. The first constraint on the 
equation of state from such an object came from the detection of gravitationally redshifted 
absorption lines in Oxygen and Iron by Cottam and collaborators [28]. By measuring a 
gravitational redshift of z = 0.35, the mass-to-radius ratio of the neutron star gets fixed at 
M/R ~ 0.15 (with M expressed in solar masses and R in kilometers). By incorporating 
additional constraints arising from Eddington and thermal fluxes, a recent analysis by Ozel 
seems to place simultaneous limits on the mass and radius of the neutron star in EXO 0748- 
676. That is, M>2.10±0.28 M and R> 13.8±1.80 km [29]. These limits are indicated 
by the black solid line in Fig. |5j An earlier determination of the spin frequency of the 
same neutron star by Villarreal and Strohmayer [30], when combined with the rotational 
broadening of surface spectral lines, yields an independent determination of the stellar radius 
of R ~ 11.5^2 5 km. This estimate, when combined with the gravitational redshift, yields 
the orange line in Fig. [5] Finally, mass-vs-radius predictions from the NL3 and FSUGold 
models are displayed in Fig. [5] The results clearly indicate the significantly harder character 
of the equation of state predicted by NL3 relative to FSUGold. This, even when both models 
predict practically identical properties for existent ground-state observables of finite nuclei. 

A critical observation that would have impacted significantly on the high-density com- 
ponent of equation of state is the one by Nice and collaborators at the Arecibo radio tele- 
scope |31j . Such (initial) observation of a neutron-star-white-dwarf binary system suggested 
a neutron-star mass of M(PSR J0751 + 1807) = 2.1±0.2 M & . This was the largest neutron- 
star mass ever reported and promised, provided that the errors could be tighten further, to 
practically pin down the high-density component of the equation of state. However, at a 
very recent conference celebrating the 40th anniversary of the discovery of pulsars in Mon- 
treal, Nice reported a significantly reduced value for the mass of PSR J0751+1807, namely, 
M(PSR J0751 + 1807) w 1.3 ±0.2 M . This revised result is denoted by the red hatched 
region in Fig. [5] and no longer invalidates any of the theoretical models under consideration. 
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FIG. 5: (color online) Constraints on the mass-vs-radius relationship of neutron stars. Displayed 
in red is the recently revised region allowed by the analysis of Nice and collaborators [31]. The 
black and orange solid lines result from the analyzes of EXO 0748-676 by Ozel [29], and Villarreal 
and Strohmayer |30j , respectively. Also shown are the theoretical predictions from the NL3 [6j [7] 
(green line) and FSUGold [TT] (blue line) models. 



IV. CONCLUSIONS 



Neutron stars are unique laboratory for the study of cold baryonic matter over an enor- 
mous range of densities. After an introduction to the "anatomy" of a neutron star, I relied 
on recent theoretical, experimental, and observational constraints to elucidate important 
features of the equation of state of neutron-rich matter. As mentioned in the Introduction, 
the only physics that spherically-symmetric neutron stars in hydrostatic equilibrium are sen- 
sitive to is the equation of state of neutron-rich matter [see Eqs. Q]. This makes neutron 
stars gold mines for the study of baryonic matter. The various constraints utilized in this 
contribution emerged from the universal behavior of dilute Fermi gases with large scattering 
lengths |25j, heavy-ion experiments that probe the high-density domain of the equation of 
state [27], and astronomical observations that place limits on masses and radii of neutron 
stars [29] [31] . On the basis of these comparisons, it was concluded that FSUGold meets all 
the challenges, even when no attempt was ever made to incorporate these constraints into 
the calibration procedure. The promise of new terrestrial laboratories (such as Facilities 
for Rare Ion Beams) together with improved observations with existent and future missions 
(such as Constellation X) offers the greatest hope for determining the equation of state of 
cold baryonic matter in the near future. 



9 



Acknowledgments 

The author is grateful to the organizers of the XXXI Symposium on Nuclear Physics for 
their kind invitation and hospitality. The author also wishes to acknowledge his many collab- 
orators that were involved in this work — especially Prof. C.J. Horowitz. This work was sup- 
ported in part by United States Department of Energy under grant DE-FD05-92ER40750. 



G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 170, 299 (1971). 

D. G. Ravenhall, C. J. Pethick, and J. R. Wilson, Phys. Rev. Lett. 50, 2066 (1983). 

M. Hashimoto, H. Seki, and M. Yamada, Prog. Theor. Phys. 71, 320 (1984). 

C. J. Horowitz, M. A. Perez-Garcia, and J. Piekarewicz, Phys. Rev. C69, 045804 (2004), 

astro-ph/0401079. 

C. J. Horowitz, M. A. Perez-Garcia, J. Carriere, D. K. Berry, and J. Piekarewicz, Phys. Rev. 
C70, 065806 (2004), astro-ph/0409296. 

G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C55, 540 (1997), nucl-th/9607039. 

G. A. Lalazissis, S. Raman, and P. Ring, At. Data Nucl. Data Tables 71, 1 (1999). 
J. D. Walecka, Annals Phys. 83, 491 (1974). 

B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

B. D. Serot and J. D. Walecka, Int. J. Mod. Phys. E6, 515 (1997), nucl-th/9701058. 

B. G. Todd-Rutel and J. Piekarewicz, Phys. Rev. Lett 95, 122501 (2005), nucl-th/0504034. 

H. Mueller and B. D. Serot, Nucl. Phys. A606, 508 (1996), nucl-th/9603037. 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 5647 (2001), astro-ph/0010227. 

C. J. Horowitz and B. D. Serot, Nucl. Phys. A368, 503 (1981). 

D. H. Youngblood, C. M. Rozsa, J. M. Moss, D. R. Brown, and J. D. Bronson, Phys. Rev. 
Lett. 39, 1188 (1977). 

J. Boguta and A. R. Bodmer, Nucl. Phys. A292, 413 (1977). 

C. J. Horowitz and J. Piekarewicz, Phys. Rev. C64, 062802 (2001), nucl-th/0108036. 

C. J. Horowitz, S. J. Pollock, P. A. Souder, and R. Michaels, Phys. Rev. C63, 025501 (2001), 

nucl-th/9912038. 



R. Michaels, P. A. Souder, and G. M. Urciuoli (2005), URL |http :77h allawe b . j lab . org/ 
|parity/prex[ 

B. A. Brown, Phys. Rev. Lett. 85, 5296 (2000). 

J. Carlson, S.-Y. Chang, V. R. Pandharipande, and K. E. Schmidt, Phys. Rev. Lett. 91, 
050401 (2003). 

G. A. Baker, Phys. Rev. C60, 054311 (1999). 

H. Heiselberg, Phys. Rev. A63, 043606 (2002), cond-mat/0002056. 

Y. Nishida and D. T. Son, Phys. Rev. Lett. 97, 050403 (2006), cond-mat/0604500. 

A. Schwenk and C. J. Pethick, Phys. Rev. Lett. 95, 160401 (2005), nucl-th/0506042. 

B. Friedman and V. R. Pandharipande, Nucl. Phys. A361, 502 (1981). 

P. Danielewicz, R. Lacey, and W. G. Lynch, Science 298, 1592 (2002), nucl-th/0208016. 
J. Cottam, F. Paerels, and M. Mendez, Nature 420, 51 (2002), astro-ph/0211126. 
F. Ozel, Nature 441, 1115 (2006). 

A. R. Villarreal and T. E. Strohmayer, Astrophys. J. 614, L121 (2004), astro-ph/0409384. 
D. J. Nice et ak, Astrophys. J. 634, 1242 (2005), astro-ph/0508050. 



10 



